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^ ABSTRACT 

The theory of the perturbative pomeron, due to Lipatov and collaborators, is 
used to compute the probability of observing parton-parton elastic scattering and 
rapidity gaps between jets in hadron collisions at very high energies. 
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The energies of the Tevatron, LHC and SSC colliders open a window on a range 
of events in a new kinematical regime of hadron scattering, the semihard regime, 
where the hadron center-of-mass energy yfs is much larger than the typical par- 
ton transverse momentum m. This complicates the experimental analysis of such 
events, since it requires jet detection in a large rapidity interval' 11 . It also makes 
harder the calculation of jet cross sections, since it involves large ratios of kine- 

[2] 

matical invariants which may appear either in the short distance cross section 

[3] 

or in the evolution of the parton structure functions , thus requiring the intro- 
duction of techniques to resum infinite classes of Feynman diagrams. Nonetheless 
the study of jet cross sections at very high energies is very interesting, since it may 
provide new insights on the dynamics of QCD processes and offer signatures of 
novel physical events' 41 . 

In order to obtain quantitative predictions of jet production in the semihard 
regime and separate it from the uncertainty involving the small x dependence 

[3] [2] 

of parton distributions , Mueller and Navelet proposed to measure the two-jet 
inclusive cross section in hadron collisions by tagging two jets at a large rapidity 
interval y and with transverse momentum of order m. These tagging jets are 
produced in a nearly forward scattering of gluons or quarks with large center-of- 
mass energy Lipatov and collaborators 1 ' 81 (BFKL) have shown that, in this 
regime, the rapidity interval y = ln(s/m 2 ) between the scattered partons is filled 
in by the radiation of additional gluons, roughly uniformly spaced in rapidity, all 
with transverse momenta of order m. The BFKL theory systematically corrects 
the lowest-order QCD result by summing the leading logarithms of s. This is done 
in three stages, as shown in Fig. 1. First, one simplifies the lowest-order QCD 
diagrams for multigluon production, shown in Fig. 1(a), for the case in which 
the emitted gluons are widely separated in rapidity. The gluon emission vertex 
is replaced by a non-local gauge-invariant effective vertex . Next, one sums the 
leading corrections to the forward amplitude with color octet exchange in the t- 
channel, as shown in Fig. 1(b). The result has the form of a Regge pole with 
an infrared-sensitive trajectory 161 . Finally, one uses this resummed, effective gluon 
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exchange to compute the elastic amplitude in the Regge limit s ^> —t with color 
singlet exchange in the t-channel, as shown in fig. 1(c). This is known as the BFKL 

ryl 

pomeron . The imaginary part of the forward amplitude is the parton-parton 
total cross section. To leading order in rapidity, the parton-parton total cross 
section and the related 2-jet inclusive cross section exhibit the energy dependence 
exp[(ap — l)y] with 



where Ca = N c = 3 is the number of colors in QCD. 

A large total cross section cr tota i for parton-parton scattering implies also a large 
elastic scattering cross section, since the total cross section is related, through the 
optical theorem, to the elastic scattering amplitude with color singlet exchange in 
the t channel. Thus both elastic scattering and jet production can be described 
by exchanging one BFKL pomeron Fig. 2(a,b). This process has been studied by 
Mueller and one of the authors' 91 . It is a higher order (a A s ) process but with energy 
dependence exp[2(o;p — l)y]. However, it leads to a final state which, at the parton 
level, contains two jets with a rapidity gap in gluon production between them 141 . 
Some fraction of these states may produce the dramatic experimental signatures of 
a large rapidity gap in secondary particle production. To understand the relation 
between rapidity gaps in hard-gluon and hadron production, we must discuss the 
potential backgrounds to these signals at the parton and hadron level. 

To analize the parton-level background, assume that we cannot detect partons 
with transverse momentum smaller than a fixed parameter /z. In this case, there 
is an additional contribution to elastic scattering from color octet exchange in the 
t channel. This is shown in fig. 1(b) where a reggeized gluon is exchanged. The 
reggeized gluon contains all the leading virtual radiative corrections and it has the 
form of a Sudakov form factor. The parameter fi fixes the scale below which soft 
gluon radiation is allowed. As ji — > 0, the contribution of the color octet exchange 
vanishes, since it is impossible to have scattering with exchange of a gluon, without 
allowing for the emission of soft gluon radiation. 
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In order to use perturbative QCD, the parameter \x must be larger than Kqcd- 
Thus we have two options: first, we can consider m ^> fi ^> Kqcd-, that is, we 
define a rapidity gap to be present if there are no jets between the tagging jets. 
We will call this case quasi elastic scattering , since it allows gluon radiation below 
the scale fx. The ratio R of the quasi-elastic to the total cross section is given by 

crtot 

where all the cross sections in (2) have been convoluted with the appropriate parton 
distributions.' 101 Alternatively, we can consider jj, = O(Aqcd)- Then at the parton 
level the color octet exchange is strongly suppressed, and only the color singlet 
exchange contributes to the cross section for producing rapidity gaps. 

At the hadron level, there is an additional consideration. Accompanying any 
hard QCD reaction in hadron-hadron collisions, there are spectator partons. These 
may produce hadrons across the rapidity interval, spoiling the rapidity gap. Thus 
in order to compute the cross section for producing a rapidity gap at the hadron 
level, we need a non-perturbative model which describes the hadron interaction 
and estimates the survival of the rapidity gap in the presence of soft spectator 
interactions' 41 . The rapidity-gap survival probability < S 2 > is defined as the 
probability that in a scattering event no other interaction occurs beside the hard 
collision of interest. This probability is most readily estimated as an average over 
the hadron-hadron impact parameter B [A] : 

n2 fd 2 Bf(B)S 2 (B) 
< S >= fdiBf(B) (3) 

where S 2 (B) is the probability that the colliding hadrons do not interact inelas- 
tically, and f(B) is the cross section for the hard collision of interest. Different 
estimates for < S 2 > in hadron collisions, based on a variety of phenomenological 
models, are presented in ref.10, where < S 2 > is estimated to be between 0.05 and 
0.2. < S* 2 > is expected to depend on the hadron-hadron center of mass energy, 
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but only weakly on the size of the rapidity gap. Then to obtain the probability of a 
scattering event with a large rapidity gap at the hadron level, we must compute the 
ratio R at /i — 0, that is, using only the singlet elastic cross section, and multiply 
it by the survival probability < S 2 >: 

R gap =<S 2 >R(fi = 0). (4) 

In this paper, we will present the asymptotic form of R(fi), and a numerical esti- 
mate for this quantity for rapidity intervals y < 15. 

For very large rapidity intervals, the one pomeron exchange approximation vi- 
olates unitarity since the elastic cross section becomes larger than the total cross 
section. Our estimate indicates that the unitarity bound for the one pomeron 
exchange approximation is at about 24 units of rapidity, when the transverse mo- 
mentum of the tagging jets m is about 20 GeV. This is far beyond the maximum 
rapidity interval available at the SSC. At lower jet transverse momenta or for larger 
rapidity intervals, R must be unitarized by the contribution from multiple pomeron 
exchange. 

Total Cross Section for Tagging Jets 

The total cross section for two tagging jets has been studied by Mueller and 
Navelet , who used the BFKL analysis of the asymptotic parton-parton total cross 
section to derive the 2-jet inclusive cross section. Following the exposition of ref. 
11, we consider the scattering of two hadrons of momenta j>a an d in the center- 
of-mass frame, with the z axis along the beam momenta, and we imagine to tag 
two jets at the extremes of the Lego plot, with the rapidity interval between them 
filled with jets. We will call these extra jets minijets. The tagging jets can be 
characterized by their transverse momenta and by their longitudinal fractions xa, 
xb with respect to their parent hadrons. It is simplest to consider the cross section 
for producing two tagging jets with transverse momenta greater than a minimum 
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value m. Then 



cIxacIxb 



(AB -> j(x A )j(x B ) + X) 



Y[ G( Xi ,m 2 ) +4/9^[Q / (x i ,m 2 ) + Q f (x il m 2 )} ■ a tot (s), 
i=A,B L / J 



(5) 



where s = 2pA -pB^AXB is the parton-parton squared center-of-mass energy, and 
&tot is the BFKL total cross section for gluon-gluon scattering within an impact 
distance of size 1/m. Eq. (5) includes the effects of quarks using the observation of 
Combridge and Maxwell that, in a process with large rapidity intervals, the leading 
contribution to any scattering process comes from gluon exchange in the crossed 
channel 1121 . The factor 4/9 in eq.(5) is the ratio of the Casimir operators Cf/Ca, 
with Cp = (N 2 — l)/2N c . The values of xa, xb should be taken to be sufficiently 
large that the parton distributions satisfy ordinary DGLAP (Dokshiter, Gribov, 
Lipatov, Altarelli and Parisi) evolution; the BFKL theory adds additional com- 
plications in the evolution of the parton structure functions when these fractions 
become small' 31 . 

The core of eq. (5) is the BFKL function at t(s), which is given by 



Ftot(y) is a dimensionless function which we will discuss below. The strong coupling 
constant is evaluated at a scale m 2 ; the running of a s is sub leading in the BFKL 
theory. Ft t is normalized to 1 at the lowest order of QCD for the gluon-gluon 
scattering cross section, integrated over the transverse momenta of the tagging jets 
with \ki±\ > m, in the limit s >> m 2 . Since Ff t is always larger than 1, it can 
be viewed through the optical theorem as the enhancement factor due to pomeron 
exchange. Ff t can be calculated from the solution of an integral equation for the 
elastic scattering amplitude, with color singlet exchange, constructed by BFKL 16 81 . 
It is given by 
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with 

u(v)=w(0,v) (8) 
the value at n = of the eigenvalue of the integral equation 

u(n, v) = ^ ty(l) - Re^(^±i + iu)} , (9) 

and ip(z) the standard logarithmic derivative of the Gamma function. The ex- 
ponential growth of Ftotiy) with the rapidity interval is associated with minijet 
production. Near u = 0, u)(v) has the expansion 



u(v) = A - Bv 2 + • • • , (10) 

with 

^ = ^ilog2, B = 1^C(8). (ID 

7T 7T 

Using (10) to expand about the saddle point at v — 0, we can see that F to t(y) 
has the asymptotic behavior 

F tot {y) ~ with z = -^j/. (12) 

V/R(3)^ 

Mueller and Navelet showed that this asymptotic form is an accurate repre- 
sentation of eq.(7) for z > 0.2. 

Elastic Cross Section for Tagging Jets 
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The high energy elastic scattering cross section with color singlet exchange in 
the t channel has been studied in ref. 9. Then it is straightforward to write down 
the elastic cross section for tagging jets, with singlet exchange, as 



d@ singlet 
dXAdXB 



{AB - j(x A )j(x B )) 

= 11 G(^,m 2 ) + (4/9) 2 ^[Q / (x,,m 2 ) + Q / (^,m 2 )] 
i=A,B<- f 



&singlet\S)- 

(13) 

The BFKL function (T s inglet(s) is the gluon-gluon elastic scattering cross section 
in the color singlet channel with the tagging jets collimated and with minimum 
jet transverse momentum m. Since two reggeized gluons are involved in the color 
singlet exchange in the t channel, in keeping into account the possibility that 
the scattering is initiated by quarks we obtain the suppression factor (Cf/Ca) 2 - 
Vsinglet(s) is given by 



CTsinglet{S) — _i) \ n J ^ m 2 r sm 9^t\V ) i 



(14) 



with 



1 singlet 



m 

167T 3 



d k A L 



f d 2 q ± d 2 q' ± f k ^(q ± ,q' ± ,y) 



(15) 



The amplitude f kA± (q±, q'±,y) represents the propagation of a BFKL pomeron 
in transverse momentum plane (fig.2(b)), with t = —k A j_. The normalization of 

/ (<Z.L)9j_>2/) * s as * n re ^- 8 so ^ at at l° west order 



f kA± (q±,q'±,y) 



S 2 (q± ~ gj) 
Q±W± ~ k A ±) 



(16) 



Lipatov 7,8 gives an integral representation of f kA± (q±, q'j_, y) in impact param- 
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eter space, for the scattering of color neutral objects. The extension to parton- 
parton scattering has been done in ref.9. The integral of f kA± (q_i_, q'^,y), over its 
transverse momenta, is' 91 



/a r 2 
d 2 q ± d 2 q' ± f k ^(q ± ,q' ± ,y) = -5- / dv " exp[a;(i/)y], (17) 
K A±_ J \ v ' 1> 

where only the leading term n — 0, for which uj(v) is given by (8) and (9), has 
been kept. As Mueller and Tang have shown, the integration over the transverse 
momenta in (17) is not infrared sensitive, i.e. the infrared divergences in qj_, which 
are present in the lowest order contribution (16), disappear when we consider the 
infinite resummation (17) of the leading logarithmic corrections. It is an interesting 
result that the pomeron trajectory has no explicit dependence on the momentum 
transfer. The dependence on t is only through the coupling constant. 

By substituting (17) into (15) we calculate F sing [ et . The integral over k_^± is 
singular and depends on the cutoff m; this gives us the same factor m~ 2 as in the 
total cross section. The enhancement factor becomes 

(2 \ 2 

/ 5^ (y2+ l)2 eXp M^J • ( 18 ) 

Fsingiet 1S normalized in order to have F s i ng i et — 1 at y — 0. Notice, though, that 
Fsingiet at y = from (18) does not correspond to gluon-gluon elastic scattering 
at the lowest order of QCD, which is instead correctly given substituting (16) into 
(15). Using the expansion of uj(v) around the saddle point v = 0, we obtain the 
asymptotics of F singlet 

/ e (41og2) 2 \ 2 

F ^ Uv) ~Ym^f>) ■ (19) 

where z is defined in (12). 
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The background to the color singlet exchange comes from the exchange of a 
reggeized gluon. This contribution is 



i=A,B 



1 1 G( Xi , m 2 ) + 4/9 Y\Qf(^ ™ 2 ) + Q f {x h m 2 )} 



f 



&octet{s,V), 



(20) 

&octet(s, A*) is the gluon-gluon elastic scattering cross section in the color octet 
channel, with the tagging jets collimated and with minimum transverse momentum 
m. We may write a oc t e t(s : //) as' 91 



&octet{s,V) = N^i^T-) 7^2 F octet(y), (21) 



c 



where F octet is the solution of the BFKL integral equation for elastic scattering 
with color octet exchange 161 



F octet = exp[2(a ff -l)j/]. (22) 
The reggeized gluon trajectory is 



a 5 = l-^ilog("), (23) 



7T \/J, 

for m 2 //! 2 ^> 1. From the expression of the gluon trajectory, we can see that 
a g — 1 < and hence F oc t e t is always less than one. As ji — > 0, or the rapidity 
interval y becomes large, & ctet(s, fx) vanishes. 

The Ratio R(fj) 

R(jJb) is the probability of having elastic scattering at the parton level, as defined 
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in (2). By substituting (20), (13) and (5) into the definition of R(fi), we obtain 

R(xA,XB,fji) = + Wf(x A )wf(x B ) asm9let (24) 

&tot crtot 

where Wf is a weight of parton distributions, 

G(x, m 2 ) + (4/9) 2 £/ [Qf(x, m 2 ) +Q f (x, m 2 )} 



Wf(x) 



G(x,m 2 ) + (4/9) EfiQA*,™ 2 ) + Q f (x,m 2 )] 



(25) 



and 

^singlet 8 / Ol s Ca \ 2 K 2 F singlet 

btot ~ N 2 -l\ 7T J 2 F tot 



' c 

&octet Foctet 



(26) 



®tot Ftot 

In the asymptotic regime, where the rapidity y is large, 

= (4 log 2) Z 



singlet 2 
~ 7T 



Ftot (iC(3)7Tz)V2 



Foctef - . e -(4 log 2+2 log(m//i)) Z 1^S)ttZ 

Ftot V 2 



(27) 



and so the octet contribution can be neglected in the asymptotic regime. The 
asymptotic expressions can be used to estimate the range of validity of the one 
pomeron exchange approximation. Taking a s ~ 0.15, which is the typical value for 
tagging jets of transverse momentum m = 20GeV, the violation of the unitarity 
bound happens at a rapidity interval y ~ 24. 

Numerical Evaluation 

We now turn to a numerical analysis of R(fi). The formulae (7), (18), (22) 
and (24) are the appropriate starting point for a numerical evaluation of the ratio 
R(n). We scale the running coupling constant from a s (m(Z)) = 0.12 using the 
1-loop evolution with 5 flavors, and take the longitudinal fractions xa, %b of the 
tagging jets to be 0.1. 
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In Fig. 3, we show the result for m = 10 — 40GeV. Even though we present 
the result from y — 0, the plots cannot be trusted in the small rapidity regime 
since the leading logarithmic approximation is no longer valid there. The dotted 
line is the contribution from color singlet exchange only, and thus corresponds to 
the choice \x = 0, while the solid line includes both singlet and octet exchanges. 
In Fig. 3(a), \x is set to lGeV. For y > 7, only the color singlet contributes since 
the octet is highly suppressed. In Fig. 3(b), /i is set to 5GeV. The suppression 
of the octet occurs now for y > 10. The two figures show the same trend in the 
asymptotic regime (y > 10) since the singlet exchange is independent of the choice 
of ii. 

The value of R(fi = 0), given by the dotted line in Fig. 3, multiplied by the 
survival propability < S 2 > (3), gives the probability of having a collision with 
a large rapidity gap (4). Bjorken M has estimated < S* 2 > to be ~ 0.1, and the 
authors of ref.10 have given estimates in the range of 0.05 to 0.2. Thus, we expect 
that, at the Tevatron, a few tenths percent of events with tagging jets will show 
rapidity gaps in hadron production. The probability of finding a gap increases with 
the rapidity interval between the tagging jets, as indicated in the figure. Though 
< S 2 > depends on the collider center of mass energy, it should depend only weakly 
on rapidity. Thus, we expect that the probability of observing a rapidity gap in 
secondary particle production between tagging jets should rise proportionately to 
the dashed curves in Fig. 3. 

Since all of the analysis above is in the leading logarithmic approximation, there 
is ambiguity in the choice of the proper scale in rapidity for which this analysis 
is valid, and so the exact value of the normalization and thus of R(fj) cannot be 
determined precisely. We need a next-to-leading order calculation, which is not 
available yet, to have a definite quantitative prediction. However, the slope of the 
curves in the asymptotic regime is free from this scale uncertainty and thus the 
experimental measurement of the ratio R(/j,)/R gap in the large rapidity-gap regime 
should give us an unambiguous determination of the survival probability < S 2 >. 
It will be interesting to see if < S 2 > depends on the width of the rapidity gap 
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between the tagging jets, when the center-of-mass energy of the colliding hadrons 
is kept fixed. 
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FIGURE CAPTIONS 



1) BFKL resummation: (a) emission of gluons; (b) construction of the octet 
amplitude; (c) construction of the singlet amplitude. 

2) BFKL Pomeron: (a) in minijet production; (b) in high energy elastic scat- 
tering. 

3) Ratio R{n) as a function of the rapidity y of the tagging jets. We choose 
xa and xb to be 0.1. The dotted line is the color singlet contribution. The 
solid line includes the contribution from both singlet and octet. We label the 
curves according to the minimum jet transverse momentum, (a) fi = lGeV; 
(b) fi = 5GeV. 



15 



